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Abstract. Optical properties of natural or designed materials are determined 
by the electromagnetic multipole moments that light can excite in the constituent 
particles. In this work we present an approach to calculate the multipole 
excitations in arbitrary arrays of nanoscatterers in a dielectric host medium. We 
introduce a simple and illustrative multipole decomposition of the electric currents 
excited in the scatterers and link this decomposition to the classical multipole 
expansion of the scattered field. In particular, we find that completely different 
multipolcs can produce identical scattered fields. The presented multipole theory 
can be used as a basis for the design and characterization of optical nanomaterials. 



1. Introduction 

The classical electromagnetic multipole expansion [Ti is a powerful tool for analyzing 
the electric and magnetic fields created by spatially localized electric charges and 
currents. Irrespective of the complexity of the charge and current distributions, the 
fields produced by them can be represented as a superposition of the fields created 
by a corresponding set of point multipoles. This correspondence provides a common 
basis for characterizing the fields radiated by localized charge and current excitations 
in arbitrary configurations. 

In optics, the multipole expansion is well suited to describe the scattering of 
optical fields by small objects. Usually, if the wavelength of the field is large 
compared to the size of the object, the scattering is described mainly by the lowest- 
order multipole, the electric dipole, while the contributions from all higher-order 
multipolcs arc considered as mere perturbations. Recently, it has been shown that in 
specifically designed optical nanomaterials [2], such as metamaterials, the contribution 
of the magnetic dipole [3] and the electric quadrupole 4 excitation to the scattering 
by the material's constituents can be made significant, which substantially affects 
the optical properties of the material and lead to extraordinary phenomena, such 
as negative refraction [5]. In certain materials, higher-order multipoles can even 
completely overshadow the electric dipole contribution [6]. Thus, it is clear that 
higher-order multipoles have to be taken into account when evaluating the macroscopic 
electromagnetic characteristics of such materials [7]. 

In order to create a material with prescribed optical properties, one should select 
the elementary unit of the material (often called the meta-atom) and optimize its 
scattering characteristics through adjusting the design. For an individual particle 
this can be done by numerically solving the Maxwell equations for the scattered field 
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and applying the multipole expansion to determine the multipoles contributing to the 
scattering [5]. However, in a material composed of a large number of such elementary 
units, each particle interacts with the fields scattered by the other particles, which 
can significantly modify the amplitudes and phases of the excited multipoles. Since at 
each point in the material the field is a superposition that contains the fields scattered 
by all the particles, the approach developed for individual particles [8] can no longer 
be used, and one should perform the multipole decomposition on the excited electric 
currents in each of the particles individually. Previously, this decomposition has been 
introduced only for a single localized electric current distribution in vacuum !lj and, 
therefore, a multipole theory suitable for analysis of nanoscatterers in an array and 
in an arbitrary homogeneous dielectric host medium was missing. In this work, we 
introduce such a theory. Our multipole expansion approach is particularly suitable 
for direct numerical implementation. 

The geometry of the scatterer determines the types of electric current modes that 
can be excited in it by light. For the description of these modes, we propose a set 
of orthogonal electric current multipoles, which are composed of elementary point 
currents in simple configurations. Each element of the resulting current multipole 
tensor reflects the strength of one of these configurations, which enables one to visualize 
the real electric current modes that will be excited in a scatterer. We complete 
our theory by deriving expressions that relate the elements of the proposed electric 
current multipole tensor to the classical multipole expansion coefficients. In particular, 
these expressions reveal i) perfectly dark multipole modes that do not create any 
electromagnetic field, and ii) electric dipole radiation produced by electric currents 
with zero net electric dipole moment. These findings provide us with additional 
freedom in the choice of the particle geometry. 

In section [2] the classical multipole expansion is adjusted to describe the 
electromagnetic field scattered by individual nanoparticles and nanoparticle arrays 
embedded in a dielectric host medium. In section [3j we map the coefficients in the 
multipole expansion to the electromagnetic fields created by sets of sub-wavelength 
current elements. Explicit mapping relations are derived up to the orders of electric 
octupole and magnetic quadrupole, which both describe third-order excitations in the 
multipole hierarchy [9 . In section |4l we summarize our results. 

2. The multipole expansion of the scattered field 

We consider a monochromatic plane wave, with the electric field amplitude Eq, angular 
frequency oj and wave vector fc, incident on a particle in an otherwise homogeneous 
lossless dielectric medium. In general, the scattered electromagnetic field can be 
written in spherical coordinates in the form of the following multipole expansion [1] 



E J2 £ i'[7r(2; + l)] 1 / 2 {-a E G,m)Vx [hP(kr)X lm (0,<l>)] 
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where Xi m and hj 1 ' are the normalized vector spherical harmonics and the spherical 
Hankel functions of the first kind, respectively. The wave number k is taken to be 
that in the surrounding dielectric which has an impedance rj. We have chosen the 
normalization of the multipole expansion such that the expressions for the scattering 
and extinction cross sections are compact. The vector functions in the multipole 
expansion form a complete basis for representing the electromagnetic field outside an 
arbitrary localized source |10j . 

Let us now assume that the scattered electric field E s is known in the region 
surrounding the scatterer (e.g., calculated numerically). Using the orthogonality 
properties of the vector spherical harmonics Xi m and the scalar spherical harmonics 
Y[ m , we can calculate the multipole coefficients from the distribution of the scattered 
electric field on any spherical surface enclosing the scatterer as 

o B (f,m) = / — / / Y? m (O,4,)r.E B (r)BmOd0d<l>, (3) 

h\ 1) (kr)E [Tr{2l + 1)1(1 + l)p/2 Jo Jo 

a M (l, m) = -~ / / Xf m (6, 0) • E s (r) sin0d0#. (4) 



/i ; (1) (fcr)£obr(2; + l)]i/2 



The above expressions for these coefficients can also be written in terms of the scattered 
magnetic field H s . In this work, we use the magnetic coefficients expressed as 

au(l,m)= \ ' 1 — / / Y l * m (9,4>)r ■ H s (r) sin 9 d9d(f>. (5) 

h\ 1) (kr)E [Tr{2l + 1)1(1 + l)] 1 /2 Jo Jo 

For scattering by a symmetric particle, all coefficients of higher order in / than some 
'max can be made equal to zero by properly selecting the origin of the coordinate 
system. The value of l max depends on both the size and the geometrical complexity 
of the particle. 

Next, we assume that the particle belongs to a large array of similar particles. 
The multipole coefficients must then be deduced from the distributions of the electric 
current density in the particles. For a single particle, such a derivation is presented in 
[1] under the assumption that the surrounding medium is vacuum. In order to allow 
for a dielectric surrounding, we define a quantity 

J s (r) = -iue [e T (r) - e lA ]E(r), (6) 

which we call the scattering current density. The electric field E = E j + J^j -^s,j 
contains both the incident field Ej (the field in the dielectric in the absence of any 
scatterers) and the field E Sy j scattered by each particle j. In ((6]), e r ,d is the real- valued 
relative electric permittivity of the dielectric and e r (r) the complex- valued relative 
electric permittivity at any coordinate r. We assume that each particle consists of 
a non-magnetic, isotropic and linear material. The assumption of discrete particles 
enables us to define a distinct scattering source current Jsj in each particle j, so that 
Js = J^j Js,y Then, starting from the ordinary macroscopic Maxwell equations, we 
derive the following equations to hold for each particle j: 

V-E Sti (r) =-|v.J 8J (r), (7) 

V-H sJ (r) -0, (8) 
V x E s ,j(r) =ikr]H S j(r), (9) 

VxH s , i (r) = ~£! s , i (r) + J s , i (r). (10) 
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In these equations, the incident field and the fields scattered by adjacent particles are 
implicitly present through Jg,j- According to (|7])-(fTU)). the introduced Jgj describes 
the effective current density that creates the scattered field of the j'th particle in the 
self-consistent solution of the Maxwell equations. 

Using ([7))- ([TO]) , we derive the scalar wave equations [1] 

(V a +* 2 )[r-J3 Bj (r)] =-iA^.J SJ (r)-i|(2 + rA)[ V . J Sii (r)], (11) 
(V 2 + fc 2 ) [r • H s j (r)] = -r • [V x J S;j (r)] . (12) 

The solutions to (fTT|) and (fT2j) are inserted into ((3]) and ((5]) to obtain the multipole 
coefficients in the form (see [1] for comparison) 

(— iV~ l kri f 

a ^ m) = E^{2i + m + i)]^ J w>*)*w 



{fc 2 r • J SJ (r) + (2 + r-)[V • J SJ (r)] }d 3 r, (13) 

OM(l,m) = gok(2 ( ^'l)|(/ + l)]V 2 fW'M'Mr • [V x J SJ (r)]d»r, (14) 

where ji are the spherical Bessel functions. While the integrations in (fTU)) and (IT4"]) arc 
over the whole space, the integrands are equal to zero everywhere outside the particle 
in question. 

The spatial derivatives of Jg,j make (IT3"|) and (|14p cumbersome, especially for 
numerical calculations. We therefore use integration by parts to cast (IT3"|) and (fT4")) in 
the form 

\l-iu2. 



,(5)0 • Jsj(r) - i7r, m (0)0 • Jsj(r)] }d 3 r, (15) 



— ; T lr< 



a M (l,m) = I ^ ^ / exp(-im(t>)ji(kr)[mi m (9)6 ■ J S j(r) 



E [tt(21 + 1)]V2 

+7i m (0)0- J SJ (r)]d 3 r, (16) 

where ^(fcr) = krji(kr) are the Riccati-Bessel functions and ^ t (kr) and ^ (fcr) are 
their first and second derivatives with respect to the argument kr. The associated 
Legendre polynomials P™ are defined as in pQ. In (fT5|) and (fT6| we have introduced 
the following functions and parameters 

21 + 1 (l -to)!] 1/2 



[1(1 + l)] 1 / 3 L 4tt (Z + to)! 

7, m (0) = ^iT (cosff), (18) 

7r im (0) = — P ; m (cos0). (19) 

Equations (IT5|) and (ITC)) yield the same multipole coefficients as © and ((H) for 
the light scattered by an isolated particle. However, equations ([3]) and ((4)) are not 
applicable to an array of scatterers. In contrast, since (fT5"]) and (fT6|) only require 
the knowledge of the total electric field inside the particles to calculate Js,j, they 
can be used to characterize the scattering of light by each particle in the array. 
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The required total electric field can be calculated, e.g., by numerically solving the 
Maxwell equations. Thus, equations (jT5j) and (jTij]) make it possible to characterize the 
optical properties of nanomaterials, in which multipoles of arbitrarily high order can 
be excited. 

For scattering of light by a single particle one can introduce the scattering cross 
section that describes the efficiency with which the particle removes energy from 
the incident plane wave into the scattered field. For our multipole coefficients the 
scattering cross section can be derived in a similar way as in [llj to become 

oo I 

C " = ££ E (2m)[K(Z,m)| 2 + |a M (Z,m)| 2 ]. (20) 

1=1 m=-l 

The terms of the series in (|2U)) allow one to determine the contribution of each 
multipole excitation to the overall scattering cross section of the particle. 

The extinction cross section can also be expressed in terms of the multipole 
coefficients. For an x-polarized incident wave of the form 

Ei(r) = xE exp(ikz), (21) 

the extinction cross section is calculated to be given by 

oo 

C C xt,x = -^5I E (2l + l)Re[ma E (l,m)+a M (l,m)]. (22) 

1=1 m=-l,+l 

In contrast to (|2T))) , the expression for the extinction cross section depends on the choice 
of the polarization and propagation direction of the incident field. For example, if the 
incident wave is y-polarized, the extinction cross section is expressed as 

oo 

, m) + mayi{l,m)\. (23) 

1=1 m=— 

The extinction and scattering cross sections can also be calculated by using the optical 
theorem [1] and calculating the total power of the scattered light. The results can then 
be compared with those obtained by using (|2"0)l . (|2"2"|) and (|2"3"|) . which can serve as an 
additional check for computations. 



3. The multipole expansion of the electric current density 

Scattering of light by a particle can effectively be seen as a process, where the 
incident light excites in the particle polarization and conduction currents that radiate. 
These currents can be decomposed into terms which we call current multipoles. To 
accomplish this decomposition, we choose to use the concept of point electric current 
elements introduced in [12 . A particle interacting with light can then be treated as a 
collection of such point elements. 

Before presenting the multipole expansion of the electric current density, we first 
consider a wire of length L that carries a time-harmonic electric current with a complex 
amplitude I. The wire is positioned at the origin of the coordinate system inside a 
homogeneous and isotropic dielectric. The oscillating current emits electromagnetic 
radiation of wavelength A into the dielectric. By assuming that L << A, we can treat 
this current-carrying wire as a point current element with the complex amplitude of 
the current density given by 

Ji(r) = IL6(r). (24) 
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This current element fully corresponds to an oscillating point electric dipole. 

Next we consider two point elements, in which the currents oscillate in opposite 
directions. The element with a current +1 is displaced from the origin in the 
positive x direction by a distance s/2, while the other element, with a current —I, is 
displaced in the opposite direction by the same amount. If s << A, we can treat this 
elementary current configuration as a second-order current element. Considering s to 
be infinitesimally small, but such that the product ILs stays finite, we obtain for the 
complex amplitude of the current density describing this second-order clement 

J 2 (r)=c(x)Ji(r), (25) 

where the operator <; is defined as 

<u) = s^- (26) 

Similar displacements can be done in the y and z directions by applying operators 
<;(y) and s{z), respectively. The point current elements of the third and higher orders 
can be obtained by sequentially applying the operator in (|26[) to the current density 
of the lowest order point element of (pM|) . In ([2"B)) . u can be chosen as x, y, z or as 
any linear combination of them. The effect of the operator on a point current element 
of a certain order can be seen as follows. The operator makes a copy of the element 
and shifts the phase of the complex amplitude of this copy by 7r radians. The original 
element is then displaced in the u direction by a distance s/2 and the copy in the 
—u direction by the same amount (see figure Q}. Finally, s is set to be infinitesimally 
small. 



(a) 



(b) 



(c) 



I 
t 



i 



Figure 1. Illustration of point current elements (a) J\ = yILS(r) , (b) 
J2 = J\ and (c) J3 = f (£)<;(y) Ji . Each red arrow represents a time-harmonic 
electric current of complex amplitude I flowing in a wire of length L. 

We describe the amplitudes of the point current elements by the following Z'th 
order current multipole moments 



(i-iy.u 



J(r) 



d 3 r 



(27) 



l— 1 terms 



where M"' is a tensor of rank I. For the orders I = 1, I = 2 and I = 3, we call 
the current multipole moments as dipole (M^ = p), quadrupole (M^ = Q), and 
octupole (7Vf( 3 ) = O) moments. For I > 2, several elements of are equal. For 

example, the current octupole moments are such that Kyz = O xzy . 
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The multipole expansion of the current density, obtained by repeatedly applying 
the operator <j to (|2~i|) , can be written in Cartesian coordinates as 

oo I— 1 l — a—1 

J(r)=k,]T E E E MW(v,a,b)v 

l — l v—x.y,z a—0 b—0 

(— lWZ— IV d Q d b d^ (a+b+1) 
' ZjxS(r), (28) 



a\b\[l - (a + b + 1)]! dx a dy b dz l - {a+b+1) 

where M^(v, a, 6), with u being equal to x, y or z, are the elements of M^> describing 
a multipole obtained by applying the operator s(u) to a ■w-oriented current element 
a times with u = x, b times with u = y and I — (a + b + 1) times with u = z. 
The coefficients in (|28|) are chosen such that the elements of the multipole tensors are 
consistent with ([27]) . 

We wish to map the elements of M"' onto the coefficients oe(Z, to) and <zm(Z, to,) 
in the multipole expansion of ([T]) and ((2]). For this we need to solve for the 
electromagnetic fields created by the current density distribution of (f2"5)) . We start by 
defining the vector potential A through 

H(r) = — V x A(r). (29) 
Mo 

In the Lorenz gauge, A satisfies the wave equation 

(W 2 + k 2 )A(r) =-/x J(r). (30) 

For a point electric current element, such as the one in (|24|) . the solution for the wave 
equation is p] 

A 1 (r) = -^h^(kr), (31) 

W 47T6 

where p = iILjuj and e = eo£ r is the electric permittivity of the surrounding dielectric. 
The multipole expansion of the vector potential is obtained by applying the operator 
^ to (f3"01 with J = J \ and A = A\, until J on the right- hand-side becomes that in 
(|28[) . Since the spatial differential operators in Cartesian coordinates commute with 
V 2 , we obtain the multipole expansion for A from (f3"U)) as 

, 3 oo l — l l — a—1 

A ^ = lLX E E E MiHv,a,b)v 

l—l v—x,y,z a—0 b—0 

i-iy-m-iy. d° d b d'-c+w) (1) 

a!6![Z-( a + fe+l)]!dx a d 2/ b dz z -( a+6+1 ) ° [ ' 

We found that the calculations can be significantly simplified by making use of a 
circular coordinate system (w,w* ,z), where w = (x + iy)/y/2 and w* is its complex 
conjugate. The unit vectors in this system are w, w* and z, where w = (x — iy)j \J2. 
These circular coordinates are linearly related to the Cartesian ones and, therefore, 
the operator <r in (|26[) still commutes with itself. Thus, the vector potential can be 
expanded as 



, 3 oo i — i L—a—L 

A ^ = lLX E EE^.M)* 

l—l i}—w^w* ,.z a—0 b—0 



(-l)'- 1 ^-!)! d a d b d'-( a+b+1 ) (1) 
a\b\[l -(a + b+ 1)]! dw a dw* b d^~( a + b+1 ) 



K l) {kr), (33) 
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E(r) 



where M ^(v,a,b) are the elements of in the circular coordinate system. The 

magnetic field corresponding to ([55)1 is given by (|2"5|) . Outside the scatterer, the electric 
field is calculated as 

iw{A(r) + ^V[V-A(r)]}, (34) 

where (|f 01) and (|3T)|) have been used. The multipole coefficients a E (Z,m) and cim(^to) 
are obtained by inserting the r components of E and H into ((31) and ([5]). For each 
multipole, the r components of E and H can be written in terms of scalar spherical 
harmonics. The orthogonality of the spherical harmonics can then be used to solve for 
a E (Z, to) and ayiil, to), without the need to manually calculate the integrals in (J3j> and 
([5]). The elements of the current multipole tensors obtained in the circular coordinate 
system can be related to the tensor elements in the Cartesian coordinate system by 
using the tensor projections. In complex coordinates, the projection of the tensor 
onto the unit vector w is ■ w* . After lengthy calculations, the following 

mapping relations are obtained (up to the electric octupole coefficients) 

a E (3, ±3) = v/l5C 3 [ ± (O xxx - O xyy - 20 yyx ) 

+ i(O yyy - O yxx - 20 xxy )] , (35) 
a E (3, ±2) = VlOC 3 [ - 20 xxz + 20 yyz - O zxx + O zyy 

± i(20 xyz + 2<3 yzx + 2<3 zxy )] , (36) 
a E (3, ±1) = C 3 [ T (30 xxx + O xyy - 40 xzz + 20 yyx - 80 zzx ) 

+ i(30 yyy + O yxx - 40 yzz + 20 xxy - 80 zzy )], (37) 
a E (3, 0) = 2V3C 3 [20 xxz + 20 yyz - 20 zzz + ZXX + O zyy ] , (38) 
o E (2, ±2) = 3C 2 [Q xx - Q yy T i(Qxy + Qy*)] , (39) 
a E (2, ±1) = 3C 2 [ T (Q xz + Q zx ) + i(Q yz + Q zy )] , (40) 
a E (2,0)= %/6C 2 [2g zz -Q xx -Q yy ], (41) 
a E (l,±l) = Ci[=FPx + ip y ] 

+ 7C 3 [ ± (O xxx + 2(9 xyy + 20 xzz — O yyx — O zzx ) 

K^yyy 2^y XX 20 yzz O xxy ^zzy)] 3 (42) 
a E (l,0)= V2C lPz 

+ 7V2C 3 [O xxz + O yyz - O zzz - 20 zxx - 20 zyy ] , (43) 
a M (2, ±2) = 7C 3 [ ± (-O xxz + O yyz + O zxx - O zyy ) 

+ i(Oxyz + Oyzx - 20 zxy )] , (44) 

<3M(2, ±1) = 7C3 [ — O xyy + O xzz + O yyx — O zzx 

T i(-O yxx + O yzz + O xxy - O zzy )] , (45) 
a M (2,0)= 7V6iC 3 [O xyz ~O yzx ], (46) 

fl M (l, ±1) = 5C 2 [ - Qxz + Ozx T i(~Qyz + Qzy)] , (47) 

a M (l,0)= 5V2iC 2 [-Qxy + Oyx], (48) 

where C x = -ifc 3 /(67re£ ), C 2 = -fc7(607re£J ) and C 3 = -ik 5 /{2107:eE a ). 

We have further verified the correctness of (f35"j) - (j48|) numerically. We 
used the computer software COMSOL Multiphysics to numerically calculate the 
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electromagnetic field created by electromagnetic point sources in subwavelength dipole, 
quadrupole and octupole configurations corresponding to each tensor element in (|35p- 
(|48p. Using the obtained electric fields in ([3]) and ((3]), we evaluated the multipole 
coefficients for each tensor element in p, Q and O. The obtained coefficients were in 
full agreement with (|35j) - (|48|) . 

The presented theory can be used as follows. For an arbitrary nanoscatterer, 
including a nanoscatterer in an array, one first numerically evaluates the fields inside 
the scatterer and uses ([6]), (fl5j) and (fT6|) to obtain the multipole coefficients aE{l,m) 
and om(^w). Then (|55|) - P5|) are used to find the essential current excitations in the 
scatterer. Each tensor element corresponds to a certain electric current mode in the 
scatterer. For example, the O xyz mode describes the current configuration obtained 
by operating with both s(y) and s(z) on an x-oriented current element. Following 
this recipe, we have, e.g., designed and characterized nanoscatterers in which incident 
light does not excite any electric dipole moment [6]. 

Let us consider some important properties of the derived current multipole 
tensors, starting with the quadrupole one. It can be seen that only 8 coefficients 
[cie(2,±2), aE(2,±l), oe(2, 0), aM(l,±l) and aM(l,0)] are used in the multipole 
expansion, whereas there are 9 elements in the Cartesian current quadrupole dyadic 
Q. This is explained by the fact that the spherically symmetric excitation with 
Qxx = Qyy = Qzz does not generate any electromagnetic field. This perfectly dark 
excitation corresponds to a radially oscillating positively charged spherical shell with 
an equal negative charge at the center. Thus, from the knowledge of the radiated 
electromagnetic field, the moments Q XXl Q yy and Q zz cannot be determined uniquely. 
However, the real currents should match the geometry of the scatterer, which allows 
one to make a unique, physically justified choice for the values of Qxx, Qyy and Q zz . 
This choice enables one to uniquely specify the excitation character. For the octupoles, 
there are 15 multipole expansion coefficients (7 electric with I = 3, 3 electric with I = 1, 
and 5 magnetic with 1 — 2), but 18 different elements in the Cartesian current octupole 
tensor. Similarly to spherically symmetric quadrupoles, the three symmetric octupole 

excitations 20xxx — ^yyx — ^zzxi ^xxy — ^^yyy — O ZZ y, and X xz — Oyyz 20 ZZ z 

are perfectly dark. One of the octupoles in each of these excitations can therefore be 
chosen arbitrarily. 




Figure 2. The electromagnetic field created by the point octupole on the left is 
indistinguishable from that created by a point dipole, when p z = 20 xxz k 2 . 
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From (|55|) - (H5|) we notice that both current octupoles and dipoles contribute to 
the same 3 dipole coefficients [aE(l,±l) and oe(1,0)] in the multipole expansion. As 
a consequence, an octupole current distribution with zero dipole moment can create 
an electromagnetic field indistinguishable from an electric dipole. For example, let 
us consider a current octupole with O — Oixxz + yyz — 2zxx — 2zyy), Q = 
and p = 0. According to (|35l) -(|48 |) . this current distribution creates exactly the same 
electromagnetic field as a dipole with a moment p — 20k 2 z. This current octupole is 
illustrated in figure [2] The equivalence between the radiation patterns of the current 
octupole and dipole can also be seen by writing the vector potential in (l3"2"j) for the 
octupole as 

A(r) = -— 20N^- + zk 2 )h { a\kr), (49) 

where we have used the fact that V 2 h^\kr) = ~k 2 h^\kr). The second term is 
precisely the vector potential of a dipole with a moment p — 20k 2 z. The first term 
is a gradient of a scalar function that does not contribute to the radiated field [see 
(|29[) ]. This finding emphasizes the importance of the relations in (|35 |) -(|48 |) when using 
multipole coefficients to describe electromagnetic excitations. The fact that spatially 
orthogonal current excitations can create the same electromagnetic fields, provides 
additional freedom to the choice of the particle geometry when designing functional 
optical nanomaterials. 

4. Conclusions 

In summary, we have introduced a theoretical approach to calculate electromagnetic 
multipole excitations in a material consisting of localized nanostructures in a dielectric 
host medium. Propagation of light through an array of such nanostructures can be 
studied numerically by solving the Maxwell equations with appropriate boundary 
conditions. Then, using our theory, the multipole excitations in the structures can 
be revealed. 

In order to obtain an intuitive picture about the real electric current excitations 
in the scatterers, we have introduced a basis of easily visualizable electric current 
multipoles such that any excitation can be represented by their linear superposition. 
The multipole expansion coefficients can be calculated numerically and used to find 
the real electric current modes in the nanostructure. The same equations can also be 
exploited in the reverse order to tailor the angular distribution and directionality of 
the scattered radiation. The theory presented in this work provides the reader with 
an exact recipe and with all necessary equations for the design and characterization 
of nanoscatterers and optical materials composed of them. 
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